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Abstract. The problem of early thermalization of matter produced in relativistic heavy- ion 
collisions at RHIC is discussed in the framework of a hybrid model that consists of the transverse- 
hydrodynamics stage followed by the standard perfect-fluid hydrodynamics stage and freeze-out. 
The two hydrodynamic regimes are connected with the help of Landau matching conditions. 
A satisfactory description of the soft hadronic observables is achieved. This indicates that the 
p5 ■ early-thermalization problem and also the HBT puzzle may be to large extent circumvented. 
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1. Introduction 

The space-time evolution of matter produced in heavy-ion collisions at RHIC can be very well 
described by the relativistic hydrodynamics of an almost perfect fluid [HE1131II1E1E1E1IE]- 
However, the use of hydrodynamics encounters the problem of early thermalization. The fast 
equilibration of matter assumed in the hydrodynamic picture is naturally understood in the 
framework of sQGP [9] or wQGP with strong plasma instabilities |1U| . On the other hand, the 
fast thermalization is difficult to explain in the framework of the microscopic models of early 
stages of collisions (e.g., in the string models). The detailed analysis of such models indicates 
that only transverse degrees of freedom undergo thermalization [TTJ . In a certain time interval 
after the collision, the system remains highly anisotropic with the transverse pressure much 
larger than the longitudinal one. The calculations using this idea (the so called transverse 
hydrodynamics) have turned out to be consistent with the data [T2 J [T3 l [T^ l [Po] . 

Despite the initial pressure asymmetries, the RHIC data suggest that the full thermalization 
eventually takes place. This observation has inspired us to construct a hybrid hydrodynamic 
model |16| which includes the step-like transition from the transverse-hydrodynamics stage to 
the perfect-fluid hydrodynamics stage. 

We assume that the system formed after the collision undergoes only transverse 
thermalization and in the proper time interval t\ < r < Ttr can be described by the tranverse 
hydrodynamics. The time t\ is the initialization time for transverse hydrodynamics, whereas 
Ttr is the transition time to the perfect-fluid hydrodynamics. During this transition the energy 
and momentum are conserved. Additionally, we require that the entropy increases. After the 
transition, for r > Tt r , the system is treated as the perfect fluid. The single freeze-out model is 
applied on the hypersurface of constant temperature T$f to calculate physical observables. The 
hydrodynamic code is coupled to the statistical hadronization Monte-Carlo model THERMINATOR 
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Figure 1. (Color online) Upper part: The comparison of the model (lines) and experimental 
(dots) transverse-momentum spectra of pions, kaons, and protons. The data are taken from Ref. 
|19| . Lower part: The comparison of the model and experimental elliptic flow coefficient V2 for 
pions+kaons and protons. The bands indicate the error of the Monte-Carlo calculations. The 
data are taken from Ref. |20|. 



|17j . which is used to calculate the soft hadronic observables: the transverse-momentum spectra 
of pions, kaons, and protons, the elliptic flow of pions+kaons and protons, and the HBT radii 
of pions. The radii are calculated with the method described in [18] . 

The idea of the initial transverse-hydrodynamics stage is used to circumvent the problem 
of very early thermalization. There exists, however, another problem connected with the 
hydrodynamic picture. The so called HBT puzzle refers to a difficulty of simultaneous correct 
description of the hadronic transverse-momentum spectra, the elliptic flow coefficient, and 
the Hanbury-Brown-Twiss interferometry data. Interestingly, the inclusion of the transverse- 
hydrodynamics stage improves the HBT results [16] . 



i 



tf 



s 

M-l 



<l) 

T) 
■H 
10 



I 

o 



(1) 

T) 
■H 
10 



n? 



6 
4 
2 

6 
4 
2 

8 
6 
4 
2 

1.3 

1.2 

1.1 

1. 



Afl nilt = 0.05 



"'out 



. STAR 

Au+Au @ Vs NN =200 GeV 

c = 20 - 30 % 

1 1 



Ar . A = 0.1 
"side 



T, = 0.25 fm 



A=0.08 
A R| = 0.02 

"long 



Ad In = 0.1 

R out/ H <;ide 



« =0.14 
n = 0.43 



T 2i = 530MeV 
T 3f = 140 MeV 



0.2 0.25 0.3 0.35 0.4 0.45 0.5 

k T [GeV] 



Figure 2. (Color online) The comparison of the model (red lines) and experimental (black dots) 
results for the pionic HBT radii. The data are taken from Ref. [21]. The model parameters the 
same as in Fig. [TJ 



2. Results 

In the upper part of Fig. [TJ we show our model results for the transverse-momentum spectra of 
pions, kaons and protons compared to the PHENIX data [19J for Au+Au collisions at the energy 
of \/snn = 200 GeV and for the centrality class 20-30%. One observes overall good agreement 
between the model predictions and the data, at the level of 16%. The model results have been 
obtained with the initial central temperature T21 = 530 MeV and the freeze-out temperature 
T%f = 140 MeV. The fraction of binary collisions in the initial profile of the density of sources 
of particle production is 0.25 and the density of transverse clusters in rapidity uq = 0.43 (our 
fitting procedure is described in more detail in [TS]). 

In the lower part of the Fig. [JJwe show our results for the elliptic flow coefficient v% compared 
with the PHENIX [20J experimental data. The pion+kaon data are reproduced at the level of 
21% while the model results for the protons overestimate the data by about 50 %. The reason 
for too large proton t>2 is the lack of hadron rescattering in the final state and/or the lack of 
viscous effects. 

In Fig. [2] we present the model HBT radii compared to the STAR [21] data. We observe very 
good agreement between the experimental data and the results taken from the model analysis. 
All radii are reproduced with the accuracy better than 10%. 



3. Conclusions 

Our analysis shows that it is possible to obtain a satisfactory description of the soft hadronic 
RHIC data in a model that does not assume a very fast thermalization of the system. Initially, 
only transverse degrees of freedom are thermalized, as suggested by the particle production 
from the fluctuating strings. Our model reproduces well the HBT radii as well as other soft 
hadronic observables. This indicates that the two main RHIC problems may be to large extent 
circumvented. The proton elliptic flow that is too large may be corrected by the inclusion of the 
shear and bulk viscosity in the standard hydrodynamic stage or by the inclusion of the hadronic 
rescattering after freeze-out. 

Recently we have also analyzed a physical scenario where the transverse hydrodynamics was 
followed by a sudden isotropization and freeze-out |14[ 115] . The present model differs from 
[HI US] by the inclusion of the extended phase described by the perfect-fluid hydrodynamics. 
The similarities between our present results and the results of Refs. [HI [15] indicate that the 
details of modeling the equilibration transition might be not very much important. However, 
the correct description of the HBT radii is achieved only if the full equilibration of the system 
is reached. 

The present approach may be generalized in several ways. In particular, the perfect-fluid 
stage may be replaced by the viscous hydrodynamics. The viscous effects may be also included 
in the transverse hydrodynamics. Naturally, the next step should be a construction of a 
continuous change from transverse hydrodynamics to perfect-fluid hydrodynamics (first steps 
in this direction have been done very recently in |22j). 

Finally we note that the large initial anisotropy of pressure may be explained only by 
introducing a very large shear viscosity, hence our model describes the initial processes which 
are out of range of application of the standard viscous hydrodynamics. 
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